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ABSTRACT 

We present continuum data from the Submillimetre Common-User Bolometer Ar- 
ray (SCUBA) on the James Clerk MaxweU Telescope (JCMT), and the Mid-Infrared 
Photometer for Spitzer (MIPS) on the Spitzer Space Telescope, at submillimetre and 
infrared wavelengths respectively. We study the Taurus molecular cloud 1 (TMCl), 
and in particular the region of the Taurus Molecular Ring (TMR). In the continuum 
data we see no real evidence for a ring, but rather we see one side of it only, appearing 
as a filament. We name the filament 'the bull's tail'. The filament is seen in emission at 
850, 450 and 160 ^m, and in absorption at 70 /im. We compare the data with archive 
data from the Infra-Red Astronomical Satellite (IRAS) at 12, 25, 60, 100 /im, in which 
the filament is also seen in absorption. We find that the emission from the filament 
consists of two components: a narrow, cold (~8 K), central core; and a broader, slightly 
warmer (~12 K), shoulder of emission. We use a radiative transfer code to model the 
filament's appearance, either in emission or absorption, simultaneously at each of the 
different wavelengths. Our best fit model uses a Plummer-like density profile and a 
homogeneous interstellar dust grain population. Unlike previous work on a similar, 
but different filament in Taurus, we require no grain coagulation to explain our data. 

Key words: stars; formation - stars: pre-main-sequence - ISM: clouds - ISM: dust, 
extinction - ISM individual: Taurus 



1 INTRODUCTION 

■ The Taurus molecular cloud is one of the nearest sta r 
forming regions, at a distance of ~140pc (|Ehasl 1 19781 ). 
As such it has often been the target of observational 
studies into the processes of molecular cloud evolu- 
tion and star formation. The large scale structure of 
the molecular cloud has been mapped using the emis- 
sion fr om CO llUngerechts fc Thaddeud [l987l : lOnishi et all 
1 19961 : [Goldsmith et al' '2005'), and also the extinction 
of starligh t (Cambrcsy 1999; Padoa n. Cambresv. fc Langerl 
l2002l : iDobashi et al|[2005l ). 

The cloud has been mapped to study 
the formation mechanism of molecular cl ouds 
jBallesteros-Paredes. Hartmann. fc Vazguez-Semadenil 
1 19991 ). and the nat ure of the turbulence t hat pervades 
them (see review bv lElmegreen fc Scald l20o3 ) . On smaller 
scales, detailed studies of the densest regions of the 
cloud have greatly increase d our knowledge of th e star- 
formation process i t self (IWard-Th ompson et alj 2000l: 
Motte fc Andrd [ioOll: lOnishi et all ,2002: Duchene et all 



Taurus differs from other nearby sites of star formation, 
in that it is forming stars in a quiescent manner. There are 
no newly formed, massive stars injecting energy into the 
cloud; and compared to more dynamic star-formation re- 
gions like Orion or p-Ophiuchus, ea ch of the star-forming 
cores are widely separated (e.g. Motte. Andre, fc Nera 



1 19981 : IWard- Thompson. Motte. fc Andrei 1 19991 ). The cloud 
also appears to have a core mass f unction which dif- 
fers from that found in other regions llOnishi et al.l I2OO2I : 
iGoodwin. Whitworth. fc Ward-Thompsonll2004l ). 

The far-infrared and submillimetre continuum 
represent an efficient means of detecting cold star- 
forming cores in molecula r clouds, even in en e rgetic 
regions such as Orion (IWard-Thompson et al.l Il995l : 



19961 ). Figure [T] 



20041 : iKirk. Ward-Thompson, fc Andrell2005h 
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iBontemps. Ward-Thompson, fc Andre 
sho ws the Taurus m olecular cloud as seen in extinction 
(Do bashi et al.|[2005l ). with the contours showing the visual 
extinctions of Av=l and Av=5. A number of the Lynds 
clouds which are associated with the Taurus molecular 
cloud are labelled. 

The majority o f the cloud has been mapped in the far- 
infrared by Spitzer (|Gudel. Padgett, fc Dougadosll2007l . for 
more details, see Section [2]2]). The region mapped by Spitzer 
is outlined in white on Figure [T] A somewhat smaller region. 



2 Nutter, Kirk, Stamatellos & Ward-Thompson 






Right ascension 

Figure 1. The Taurus molecular cloud seen in extinction (jPobashi et al.l[2005^ . The Av=l and Av=5 contours are shown. TMCl is 
marked, as are a number of Lynds dark clouds llLyndsHlQeaT The outlines of the regions mapped by MIPS at 160/xm (large box), and 
SCUBA at 850/im (small box), are shown. The SCUBA data are centred around TMCl. 



centred on Taurus molecular cloud 1 (TMCl), has been 
mapped by SCUBA on the JCMT (see Section [23 . This 
smaller region is also outlined in white on Figure [1] The 
densest part of t he cloud is T MCl. Together with Heiles 
Cloud 2 (HCL2; iHeilesI 119681) this makes u p the Taurus 
Molecular Ring (TMR) (see lToth et~ai1l2004l '). 

Figures [21 a) and (b) show the region around TMCl at 
160 and 850 /im respectively, with the brighter components 
of Heiles Cloud 2 and TMCl marked. In addition, the ap- 
proximate locus of the TMR is also marked. For more details 
about these data, see Section [2] However, we note here that 
the molecular ring is not really seen in the continuum. In- 
stead, we see a filament in the eastern half and a series of 
point-like and extended sources scattered roughly around 
the western half. We therefore challenge the existence of the 
ring as a single coherent entity, as that is not what is seen 
in these data. We here name the filament 'the bull's tail'. 

In this paper, we present a study of a portion of 
TMCl using infrared and submillimetre data between 12 
and 850 fim. This range of wavelengths straddles the peak 
of the spectral energy distribution (SED) for the cold dust 
that is found in the densest regions of molecular clouds. In- 
frared and sub-mm data are therefore a powerful tool for 



breaking the degeneracy of temperature and density that is 
present when only sub-mm data are considered. 

2 OBSERVATIONS 
2.1 SCUBA Data 

The submillimetre data presented in this study were ob- 
tained using the Submillimetre Com mon User Bolometer 
Array (SCUBA - [Holland et~al]|l999l ') on the James Clerk 
Maxwell Telescope (JCMT). SCUBA takes observations at 
450 and 850 j^m simultaneously through the use of a dichroic 
beam-splitter. The telescope has a resolution of 8 arcsec at 
450 fim and 14 arcsec at 850 fim. The data presented here 
were acquired from the JCMT data archive, operated by the 
Canadian Astronomy Data Ce ntre. A sub-set of th ese data 
has been published previously jSchnee et al.ll20o"7h . though 
all of the data have been re-reduced by us using a consistent 
method for the purposes of this study. 

The observations were carried out over 12 separate 
nights between December 1998 and January 2005 using the 
scan-map observing mode. A scan-map is made by scanning 
the array across the sky, using a scan direction of 15.5° from 
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Figure 2. A close-up of the TMCl region at (a) 160 /im, and (b) 850 /im. The brigli t components of TMCl and HCL 2 are marked 
(TMCl-CP and TMC1-NH3 refer to tlie cyclopolyyne and ammonia peaks respectively - lOlano. Walmslev. fc WilsonlfigSSh . Some bright 
IRAS sources are also marked. The approximate locus of the TMR is shown as a dashed line. The ring structure is not very clear in the 
continuum images. Instead we see a filament coincident with the eastern portion of the ring, and simply a collection of point sources and 
clumps in the western half. We here name the filament 'the bull's tail'. 



the axis of the array in order to achieve Nyquist sampling. 
The array is rastered across the sky to build up a map sev- 
eral arcminutes in extent. 

Time-dependent variations in the sky emission were re- 
moved by chopping the secondary mirror at 7.8 Hz. The 
size of a scan-map is larger than the chop throw, therefore 
each source in the map appears as a positive and a negative 
feature. In order to remove this dual-beam function, each re- 
gion is mapped six times, us ing chop throw s of 30, 44 and 68 
arcsec in both RA and Dec jEmersonlllOoH ). The dual-beam 
function is removed from each map in Fourier space by di- 
viding each map by the Fourier transform of the dual-beam 
function, which is a sinusoid. The multiple chop-throws al- 
low for cleaner removal of the dual beam function in Fourier 
space. The maps are then combined, weighting each map 
to minimise the noise introduced at the spatial frequencies 
that correspond to zeroes in the sinusoids. Finally the map 
is inverse Fourier transformed , at which point it no lon ger 
contains the negative sources l|jenness fc Lightfootllioool l. 

The submillimetre zenith opacity at 450 and 850 /im 
was determined using the 'skydip' method and by com- 
parison with polynomial fits to the 1.3 mm sky opacity 
data, measured at the Caltech Submillimeter Observatory 
I Archibald et al.l 120021 '). The sky opacity at 850 /im varied 
from 0.18 to 0.59, with a median value of 0.27. These cor- 
respond to a 450 /im opacity range of 0.81 to 3.5, and a 
median value of 1.4. 

The data were reduced in the normal way us- 
ing the SCUBA User Reduction Facility, SURF 
jjenness fc LightfootI |2000| ). Noisy bolometers were re- 
moved by eye, and the baselines, caused by chopping onto 



sky with a different level of emission, were removed using 
the MEDIAN filter. Calibration was performed using obser- 
vations of the planets Uranus and Mars, and the secondary 
calibrator CRL618 (SandcU 1994) taken during each shift. 
We estimate that the absolute calibration uncertainty 
is ±5% at 850 fj,m and ±15% at 450 fim, based on the 
consistency and reproducibility of the calibration from map 
to map. The resolution of the SCUBA data was degraded 
using a Gaussian smoothing routine to a 45 arcsec beam, 
to allow a direct comparison of these data with the Spitzer 
data. 



2.2 MIPS Data 

The Multiband Imaging Photometer for Spitzer (MIPS; 
iRieke et al.ll2004^ is the far-IR camera on board the Spitzer 
Space Telescope jWerner et al. 20041. It operates at 24, 70 
and 160 /im, and has a diffraction-limited resolution of 6, 18 
and 40 arcsec respectively at the three wavelengths. 

Spitzer legacy surveys have been used to map most of 
the star-forming regions within 0.5 kpc and hence produce a 
complete mid-infrared record of the nearby large molecular 
clouds. These legacy surveys include 'Cores to Disks' (c2d 
- lEvans et"al] 120031. the 'Spitzer Gould Belt Legacy Sur- 
vey' (lAUen et a l. 2007) , and the 'Taurus-2 S pitzer Legacy 
ProiectlH (ICMel. Padgett, fc Dougadosll2007h . The Spitzer 
archive data used here were observed as part of the 'Taurus- 
2 Spitzer Legacy Project'. 



^ |http : //spider ■ Ipac . caltech. edu/staff/dlp/taurus/| 
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Figure 3. An even closer view of the bull's tail filament at (a) 850 /im, (b) 450 fim, (c) 160 fim, and (d) 70 ^m. Tfic filament is seen 
in emission at the tiiree long wavelengths, and in absorption at 70 fira. Each of the maps has been smoothed to match the 45-arcsec 
angular resolution of the MIPS 160-/im data. The remaining differences in the width of the filament are therefore intrinsic to the source 
itself. The filament is narrowest at 850 and 70 /um, where there is good agreement between the emission at 850 and absorption at 70 /im. 
The filament appears broader at the other two wavelengths. 
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Figure 4. One-dimensional cuts through the data perpendicular to the filament at each of the eight wavelengths studied. Each cut has 
been normalised to its own peak. 



The LEOPARD software tool (|LE0PARD| l2005h was 

used to search the Spitzer data archive for MIPS scan- 
map data coincident with the position of the SCUBA 
data. 160-/^m unfiltered basic calibrated datasets 1122- 
6112/6368/6624/6880 from pipeline version S14.4 were iden- 
tified and downloaded. The data were filtered to remove 
the frames taken directly after stimulator-flashes (calibra- 
tion events that can leave bright residuals). The data 
were re-gridded with the MOPEX software package into 
an im age with 15-arcsec-diameter pixels (Makovoz &: Khan 
l2005l ). This image was then iteratively cleaned and 
filled to reject spurious pixel a rtefac ts as described in 
iKirk. Ward-Thompson, fc Andrei (|2007| ). 

Rather than process the 70- /im data ourselves we have 
made use of data released by this Legacy Project. These 
data come from the same original dataset, but have been 
processed to an enhanced level by the Legacy Project team. 
The data were originally reduced with pipeline version SI 1.1, 
and finally smoothed to the same 40-arcsec resolution as the 
160-/^m data. 

The chopping that is employed by the JCMT to remove 
the effects of the atmosphere also has the effect of limit- 



ing the sensitivity to spatial scales larger than a few times 
the largest chop throw. Spitzer does not need to chop, and 
so does not suffer from this. In order to directly compare 
the MIPS and SCUBA data, synthetic chops with the same 
properties as the SCUBA chops were added to the MIPS 
data using the SURF routine ADD_DBM. These chops were 
then removed using the same method as that employed for 
the SCUBA data. The SCUBA and MIPS data therefore 
have the same angular resolution, and are sensitive to struc- 
tures on the same spatial scales and can thus be compared 
directly. 

2.3 IRAS Data 

The Infrared Astronomical Satellite (IRAS: IWheelock et al.l 
Il994h was a 60-cm space-borne telescope which carried out 
all-sky surveys at 12, 25, 60 and 100 /^m. IRAS had an 
angular resolution of 20, 20, 60 and 120 arcsec respec- 
tively at the four wavelengths. ISSA data for TMCl were 
retri eved from the IRAS archive using the skyview inter- 
face (|McGlvnn. ScoUick. fc WhitdflOOj ) . In order to remove 
large scale gradients which were present in the IRAS data. 
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they were smoothed using a 5-arcmin Gaussian function. 
This smoothed map was subsequently subtracted from the 
data, thus removing the large-scale structure, and also mak- 
ing the data directly comparable with the chopped data. 



3 RESULTS 

Figures [3l;a)-(d) show the bull's tail at 850, 450, 160 and 
70 jj,m respectively. The three long wavelength maps show 
the filament in emission, whereas the 70 jj,m map shows the 
filament in absorption against the warmer background of the 
surrounding cloud. This is therefore clearly a dense struc- 
ture, for it to be seen in absorption at a wavelength as long 
as 70 /xm. 

Figure |3] shows that the filament is much narrower at 
850 than at 450 or 160 /im. In addition, the absorption at 
70 fim is also very narrow, and matches quite well to the 
850-/xm emission. We interpret this as the filament having a 
cold dense inner region, where the extinction is the highest, 
surrounded by a warmer, less dense outer 'jacket'. This is 
explored further in Section |4TT] 

Figure |4] shows a series of one-dimensional cuts through 
the data of Figure[3l in a direction perpendicular to the long 
axis of the filament. In addition. Figure |3] shows the equiv- 
alent cuts through the IRAS data at 12, 25, 60 & 100 /im. 
Figure |4] shows that there is a transition at approximately 
100 nm, longward of which the filament is seen in emission. 
At wavelengths shorter than 100 fim the filament is seen in 
absorption. At 100 fim there is a hint that the filament is 
seen in absorption, although the data are quite noisy. Fig- 
ures |4ja)-(c) also confirm that the filament is narrower at 
850 fim than at 450 or 160 fim. We investigate this effect in 
Section im 



4 ANALYSIS 

4.1 Empirical Model 

Figure [5] shows the 850- and the 160-/im profiles overlaid. 
The two profiles have been normalised to a blackbody with 
a temperature of 12 K. To do this, the observed 850- and 160- 
Hm profiles were divided by the flux-density which would be 
emitted by a 12 K blackbody with a 160-/im flux-density 
equal to the observed value at the filament centre. The two 
profiles would have the same peak value of 1 if the data were 
due to a single emitting source at 12 K. We use a blackbody 
function in this instance, as it has fewer free parameters than 
the modified blackbody used in Section [4.21 

In the outer parts of the filament, there is quite a good 
agreement between the two profiles, indicating that the same 
dust is emitting at both wavelengths, and that the 850-/im 
data also match this model quite well in the 'shoulders' of 
the profile. However, there is an excess of emission at 850 /im 
in the centre of the profile. This shows that there is a second 
dust component in the centre of the filament which is colder 
than 12 K, which therefore emits more strongly at 850 fim. 
In Section [S] we discuss specific models of the filament, but 
first we attempt a simple spectral energy distribution (SED) 
analysis of the two components. 
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Figure 5. The 850- and 160-/xm cuts normalised to a 12-K black- 
body (see text for details). Note that the two curves match up 
well in the 'shoulders' of the two profiles, but that the 850-/xm 
data have an additional narrow component of extra emission on 
the peak of the filament. 

4.2 Temperature Analysis 

We measured the total flux density from the filament at each 
of the three wavelengths in which it is seen in emission. We 
also separated out the flux density of the central component 
from that of the outer shoulder, or jacket, by fitting two 
components to the various profiles, such as those seen in 
Figure [5] 

Figure [6j a) shows the flux densities of the filament cen- 
tre (solid squares) and the outer shoulder (hollow squares) 
together with two modifie d blackbody fits to the data in the 
norm al manner ( see, e.g . 'Ward-Thompson, Andre, & Kirl^ 
l2002h . Following iKirk. W ard-Thompson. & Andre (200!$^ 
the parameters j3 and the critical wavelength (Ac) are con- 
strained to be equal to 2 and 50 ^m respectively. The best 
fits to the data are found to be 12 K for the outer filament, 
and 9 K for the filament centre. 

We assume when we look along a line of sight towards 
the centre of the filament, that as well as the cold central fil- 
ament, we also see the outer component (in front of, and be- 
hind the central filament). Therefore, to model the emission 
from the central component, we must also take account of 
the outer component. Hence we must fit a two-temperature 
SED to the emission from the central component. 

Figure IHl^b) shows the flux densities of the central com- 
ponent of the filament fitted by two modified blackbodies. 
We constrain the warmer component to be 12 K, as dis- 
covered in Figure Oa), and fit the residual emission with 
a second component. In this way we find that the best fit 
temperature for the cold dust is actually 8 K. The resulting 
two-part SED is plotted on Figure |S][b). Hence we see that 
the best fit to the data from a simple two-component SED 
analysis is for a central cold filament at 8 K, surrounded by 
a broader and warmer component at 12 K. 

4.3 Comparison with previous work 

iToth et all l|2004l ') used the IRAS 60- and 100-^im data, to- 
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Figure 6. (a) Flux densities of the central component of the 
filament (filled squares) and the outer component (open squares) 
at each of the three wavelengths at which the filament is seen in 
emission. Spectral energy distribution fits in the form of modified 
blackbody curves are also shown for both the central filament 
(solid line - 9 K) and the outer filament (dashed line - 12 K). (b) 
Flux densities for the central component of the filament only. A 
two-component SED for the central part of the filament is also 
shown. The dashed line represents the same 12-K fit as before, 
but now the dotted line shows a fit to the residual flux densities at 
the two longer wavelengths, which corresponds to an 8-K modified 
blackbody. The solid line shows the sum of the two. 



gether with 200-/im data from the Infrared Space Obser- 
vatory (ISO) to determine the temperature of TMCl-CP, 
which is roughly coincident with the filament that we are 
studying (see Figure [2]). 

They estimated a temperature of 11.4 K, which is 
consistent with the 12 K that we have me asured in the 
outer part of the filament. iToth et al.l l|2004l ') also found a 
similar tem perature from studies of NH3(1, 1) and (2,2). 
IStepnik et all (|20o3 ') used the baUoon-borne PRONAOS 
telescope to observe the LI 506 cloud in Taurus (see Fig- 
ure[ T]). They mea s ured a temperature of 12.1 K. 

Schnee et all ((2003) observed the TMCl-C region us- 



ing SCUBA data at 450 and 850 fim together with 1.2-mm 
data from the Institute de Radio Astronomie Millimetrique 
(IRAM) 30-m telescope. They found the temperature to be 
in the range of 6 - 16 K. All of these independent estimates 
are consistent with our findings. 



4.4 Column density 

The measured values of the 850 /im flux density at the centre 
of the filament and in its outer part are 0.082 Jy/beam and 
0.022 Jy/beam respectively. The fiux density of the outer 
part was measured 100 arcsec from the centre of the fil- 
ament. From these measurements, together with the tem- 
peratures determined above, we can calculate the column 
density of material along the two lines of sight. 
The mass per beam can be calculated using: 



K850-B850 



(1) 



where Sgso is the 850 nm flux density, D is the distance to 
the fllament, K850 is the mass opacity of the gas and dust, 
and Bs5o,T is the value of the blackbody function at a wave- 
length of 850 at t emp erature T. We assume a distance 
to TMCl of 140 pc mias 197 i), and a mass opacity of 
0.01 cm^g~^ (see [ Andre. Ward- Thompson, fc Barson- 



19931; 



Andre. Ward- Thompson, fc MottJ [Tggi 



jsonyl 
19961 : 



Ward- Th ompson. Motte. fc Andrei 1 19991 . for a detailed 
discussion of both this value of Kgso in particular and of 
this method of obtaining masses in general). 

The column density, N{H2), is then calculated using: 



N{H2 



Mb, 



(2) 



where rbeam is the radius of the beam at the distance of 
TMCl. This yields column densities of 3.8 x 10^^ cm"^ and 
4.4 X lO^'^ cm~'^ for the central part and the outer part of 
the filament respectively. Hence we see that the centre of the 
filament has roughly an order of magnitude higher column 
density than the outer part. 

We can also calculate the column density using the 70- 
^m data, where the filament is seen in absorption, using 
the method described by Bacmann et al. (2000 - see also 
Ward- Thompson et al. 2006). This compares the measured 
intensity where the absorption is the strongest (/on), with 
the intensity off the filament (loff)- A third measure {I fore) 
away from the region is also required to account for the 
zodiacal hght and other foreground emission between us and 
the molecular cloud, as well as the large-scale background 
Galactic emission. The optical depth (ta) of the fila ment is 
given by the following (|Ward-Thompson et al.|[2006l ): 



e ^ — . 

loff — I fore 

The column density is calculated using: 



(3) 



(4) 



where o"a is taken to be 2.6 x 10"^'' cm^ at 70 fim 
JPraine fc Led [l984l . see their figure 9). Measured values 
of Ion, loff and I fore are 50.91, 51.13 and 49.51 MJysr"^ 
respectively, which lead to an estimate of the column den- 
sity of 5.6 X 10^^ cm~^. This is very close to the value of 
3.8 X 10^^ cm~^ calculated above from the 850-/im data in 
emission. The fact that two different data-sets from two dif- 
ferent telescopes, and two different methods of calculation, 
give essentially the same answer give added confidence to 
our estimates. 
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5 RADIATIVE TRANSFER MODELLING 

In this section we employ more sophisticated techniques to 
model all of the data, in order to understand the density 
structure of the filament in more detail. We use a radiative 
transfer code (Stamatellos & Whitworth 2003) to model the 
data, using a variety of parameter values to obtain the best 
fit simultaneously at all wavelengths. In each case, the fila- 
ment is approximated by a cylinder, which is heated exter- 
nally by the interstellar radiation field (ISRF). To determine 
the properties of the filament a large number of radiative 
transfer simulations were performed, and a best-fit solution 
was found. 

We assume that the dust consists of carbona- 
ceous and silicate g r ains with size distribution from 
IWeingartner fc Draind (|200ll ). These propert ies are appro - 
priate for the dust in the ISM {Rv = 3.1; see lDraineir2003h . 



5.1 Monte Carlo radiative transfer 

The radiative transfer calculations were performed using 
PHAETH ON, a 3D Monte Carlo radiative t ransfer code de- 
veloped bv lStamatellos fc WhitworthI l|2003l ). The code uses 
a large number of monochromatic luminosity packets {L- 
packets) to represent the radiation sources in the system. 
The L-packets are injected into the cloud and interact (are 
absorbed, re-emitted, scattered) stochastically with it. 

If an L-packet is absorbed its energy is added to the 
local region and raises the local temperature. To ensure ra- 
diative equilibrium the L-packet is re-emitted immediately 
with a new frequency chosen from the difference between 
the local ce ll emissivity before and after the absorpti on of 
the packet jBiorkman fc Woodll200ll:lBaes et al.ll2005h . For 
more details, see IStamatellos fc WhitworthI (|2003l V 

We modelled the filament using a symmetric cylinder, 
and hence the code used here is optimised for the study 
of systems with axial symmetry. The filament is typically 
divided into ~ 4000 cells. The number of cells used is chosen 
such that the density and temperature differences between 
adjacent cells are small. 



5.2 Heating of the filament 

Our simple modelling above has shown us that the centre of 
the filament is colder than the outer parts. Thus we deduce 
that there are no heating sources embedded within the fil- 
ament, and hence the filament is starless. In order to verify 
this, Spitzer MIPS 24 ^m and Spitzer IRAC dat a for the fila- 
ment were searched for YSOs and protqstars (c.f. lAllen et al.l 
l2004l : iHarvev et allbOOTl : iRebuU et al.]|2007l ). With the ex- 
ception of a single class I source coincident with the out- 
skirts of the filament, there are no observed internal heating 
sources. The only heating source considered is therefore the 
interstellar radiation field (ISRF) . For this ra diation field 
we adopt a revised version of the iBlackl (|1994| ) interstellar 
radiation field (BISRF). The BISRF consists of radiation 
from giant stars and dwarfs, thermal emission from dust 
grains, cosmic background radiation, and mid-i nfrared emis- 
sion f rom transiently heated small PAH grains (lAndre et al.l 
l2003l ). 

This radiation field is modulated by the ambient cloud 
around the filament. Hence the incident radiation field on 
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Figure 7. A Plummer-like analytic density profile, with 
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the filament is enhanced at FIR and longer wavelengths, 
and attenuated at shorter wavelengths. The L-packets rep- 
resenting the ambient radiation field (typically a few 10^" 
packets) are injected from the outside of the virtual ambi- 
ent cloud surrounding the filament with injection points and 
injection directions chosen to mimic an isotropic radiation 
field. 



5.3 Plummer-like geometry 

The geometry we assume is a c ylinder with a Plummer-like 
density profile (|Plummerlll91ll) . So we use axial symmetry 
with the density being approximately fiat in the centre of 
the filament and dropping as a power-law in the envelope of 
the filament, according to the Plummer relation: 



Rflat 



(5) 



where r is the distance from the axis of symmetry of the 
cylinder, n{H2)fiat is the density at the axis of symmetry, 
Rflat is the extent of the region in which the density is 
approximately uniform, and rj is the power-law slope at large 
values of r. This analytic form is illustrated in Figure [7] The 
filament is surrounded by a virtual ambient cloud, which 
has a uniform density n{H2)amb- The role of this cloud is to 
modify the ambient radiation field that heats the filament 
externally. 

We use the term Plummer-Uke because a true Plummer 
profile has = 5, whereas we treat this as a free parame- 
ter. A Plummer-like profile was selected because it provides 
a high-density inner region, which decreases to large radii, 
with a small number of parameters. The model can produce 
images of the source at any requested wavelength, convolved 
to any resolution. From these images we can produce one- 
dimensional cuts perpendicular to the cylinder for compar- 
ison with_the_data; 

IWhitworth fc Ward-ThompsonI (1200 ll ) used a 
Plummer-like profile to model prestellar cores as they 
collapse to form protostars. They found that this model 
successfully predicted both the observed infall rates, and 
the relative lifetimes of the different evolutionary stages. 
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Figure 8. The cuts through the data from Figure |3] (bold lines), with overlaid (thin lines) the best fit profiles from the radiative transfer 
model. 



The following parameter space was explored. Rfiat = 
5 X 10^ - 1.5 X 10* AU, n[H2)fiat = 7 x 10'' - 2 x 10^ 
cm~'^, = 1 — 4. In addition, the filament was embedded 
in a diffuse medium, with a visual extinction varying from 
Av =0.01-0.55. 

Figure |S] shows the best-fit results for this model. The 
cuts through the model filament are compared to those 
through the data (see Figure [l]). The parameters for this 
model are Rfiat = 10* AU, n{H2)fiat = 1-75 x 10^ cm"^, 
rj — 3 and Av in the external cloud equal to 0.39. The IRAS 
data are not included in this figure as they are of insufficient 
quality for detailed modelling at this scale. 

This model shows a good fit to the filament in emission 
at 850 and 160 fim, and fits one side of the asymmetric 
profile of the filament at 450 ^m. The model also shows the 
filament in absorption at 70 fim. While this is not necessarily 
a unique fit to the data, it shows that the filament can be 
modelled with a fairly simple geometry and uniform grain 
properties throughout. 



This is in contrast to the result of lStepnik et al.l (|2003l ) , 
who need to invoke a non-uniform grain distribution in order 
to fit the far-infrared de nsity profiles of a similar filament 
in Taurus. IStepnik et al.l conclude that grain aggregation is 
important in the central 3 x 10* AU, and is required to model 
the emission. We can explain the emission from our filament 
without a need for a change in grain properties. 

The radiative transfer model allows us to measure the 
temperature of the dust grains as a function of position. This 
temperature does not vary with grain size, but is the same 
for all grains at a specific position. The temperature profile 
of the best-fit model is shown in Figure (9] This gives us an 
independent temperature estimate for comparison with the 
SED fitting carried out in Section [4.21 The model predicts 
a temperature of 12 K for the outer component of the fila- 
ment, in agreement with the temperature determined from 
the SED fits. The central temperature is slightly higher than 
the best fit SEDs (9.8 K compared to 8 - 9 K). 
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Figure 9. The temperature profile from the best-fit model. The 
vertical dotted line shows the radial distance of the outer compo- 
nent of the filament (see section |4]2]) 



6 CONCLUSIONS 

We have presented submillimetre data from SCUBA on the 
JCMT and far-infrared data taken by Spitzer and IRAS, for 
a filament in TMCl which we have named the bull's tail. 
TMCl is the densest part of the Taurus molecular cloud 
region. The filament is seen in emission at all wavelengths 
longer than 100 fim, and in absorption at all shorter wave- 
lengths. 

The filament is significantly narrower at 850 ^m, which 
we attribute to a narrow, dense inner filament, running 
through the centre of the broader and warmer main fila- 
ment. We have shown intensity cuts through the filament 
at wavelengths of 850, 450, 160, 100, 70, 60, 25 and 12 fim, 
which are consistent with this scenario. 

We have fitted SEDs to both the inner and outer com- 
ponents of the filament and find that the temperature de- 
creases from 12 K in the outer part to 8 K at the centre. We 
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then used a radiative transfer code to fit the profile of the 
filament. Wc modelled the filament using a Plummor-like 
density profile with constant grain parameters throughout, 
and find good fits to all wavelengths. 
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